Further matters in space-time geometry: / (R, T, R UV T^ U ) gravity 
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We consider a gravitational theory in which matter is non-minimally coupled to geometry, with 
the effective Lagrangian of the gravitational field being given by an arbitrary function of the Ricci 
scalar, the trace of the matter energy-momentum tensor, and the contraction of the Ricci tensor with 
the matter energy-momentum tensor. The field equations of the theory are obtained in the metric 
formalism, and the equation of motion of a massive test particle is derived. In this type of theory the 
matter energy-momentum tensor is generally not conserved, and this non-conservation determines 
the appearance of an extra-force acting on the particles in motion in the gravitational field. It is 
interesting to note that in the present gravitational theory, the extra-force explicitly depends on 
the Ricci tensor, which entails a relevant deviation from the geodesic motion, especially for strong 
gravitational fields, thus rendering the possibility of a spacetime curvature enhancement by the 
RyujT^" coupling. The Newtonian limit of the theory is also considered, and an explicit expression 
for the extra-acceleration which depends on the matter density is obtained in the small velocity 
limit for dust particles. We also analyze in detail the so-called Dolgov-Kawasaki instability, and 
obtain the stability conditions of the theory with respect to local perturbations. A particular class of 
gravitational field equations can be obtained by imposing the conservation of the energy-momentum 
tensor. We derive the corresponding field equations for the conservative case by using a Lagrange 
multiplier method, from a gravitational action that explicitly contains an independent parameter 
multiplying the divergence of the energy-momentum tensor. The cosmological implications of the 
theory are investigated in detail for both the conservative and non-conservative cases, and several 
classes of exact analytical and approximate solutions are obtained. 

PACS numbers: 04.30.-w,04.50.Kd,04.70.Bw 



I. INTRODUCTION 

The recently released Planck satellite data of the 
2.7 degree Cosmic Microwave Background (CMB) full 
sky survey [l|, l2| have generally confirmed the standard 
ACDM (ACold Dark Matter) cosmological paradigm. A 
major goal of the Planck experiment was to test the 
ACDM model to high precision, and identify possible ar- 
eas of tension. In fact, an interesting result from the fits 
of the basic ACDM model to the Planck power spectra is 
the lower than expected value of the Hubble constant Hq, 
Hq = 67.3 ± 1.2 km/s/ Mpc. The Hubble constant can 
be tightly constrained by CMB data alone in the ACDM 
model [2| . The Planck data has also further constrained 
the parameters of dark energy, a possible cause of the 
late-time cosmic acceleration. Indeed, a central problem 
in present day physics is to elucidate the nature of dark 
energy, thought to be driving the accelerated expansion 
of the Universe. Perhaps the most straightforward ex- 
planation for dark energy is the presence of a cosmolog- 
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ical constant. An alternative is dynamical dark energy 
|3,|4J], usually assumed to be a very light scalar field, hav- 
ing a canonical kinetic energy term, and being minimally 
coupled to gravity. The cosmological constant A has an 
equation of state w = pj p = — 1, where p and p are the 
effective thermodynamic pressure and energy density as- 
sociated with A, while scalar field theories usually have 
time varying equations of state with w > — 1 [2| . 

The CMB alone does not strongly constrain the dark 
energy equation of state parameter w, due to the two- 
dimensional geometric degeneracy present in dark en- 
ergy models. However, this degeneracy can be broken 
by combining the CMB data with lower redshift dis- 
tance measurements [2j. By combining the Planck data 
with the measurements of Hq in yj], the authors pro- 
vide an equation of state parameter of dark energy given 
by w = —1.24^0^, which is off by more than the 2cr 
compared to id = —1 [S|. The Planck data combined 
with the Baryon Acoustic Oscillations (BAO) data give 
w = — 1.131q'25 0- Therefore, presently there is no con- 
vincing observational evidence that could clearly estab- 
lish the nature of dark energy. Moreover, the accelerated 
expansion of the Universe (see [5[ and references therein) , 
the virial mass discrepancy at the galactic cluster level 
and the galaxy rotation curves [6| as well as other cos- 
mological observations suggest that the standard gen- 
eral relativistic gravitational field equations, based on 



the Einstein-Hilbcrt action S — J (R/2 + L m ) ^J~gd A x, 
where R is the scalar curvature, and L m is the matter La- 
grangian density, cannot describe the Universe at large 
scales, beside passing the Solar System tests. From a 
cosmological viewpoint, this amounts to introducing, by 
hand, the dark matter and dark energy components in 
the theory, in addition to ordinary matter and energy. 

Another possibility is to modify the basic structure of 
the Einstein-Hilbcrt action in the hope that such a modi- 
fication could naturally explain dark matter and dark en- 
ergy, without resorting to some exotic forms of matter. 
Initially, the interest in the extension of the Einstein- 
Hilbert action was focused on the modification of the 
geometric part of the action. One of the interesting re- 
search avenues is the introduction of higher order terms 
to the gravitational field action through the substitution 
of the Ricci scalar with a generic function f(R). The so- 
called f(R) gravitational theories were first introduced 
in [Tj, and later used to find a non-singular isotropic de 
Sitter type cosmological solution [8| ■ For detailed reviews 
of f(R) theories see, for instance, [3|. The most serious 
difficulty of f(R) theories is that in general, these theo- 
ries seem incapable of passing the standard Solar System 
tests Pi|- However, there exists some theories that can 
accommodate this problem |ll| . On the other hand, us- 
ing phase space analysis of the specific model involved, 
it has been shown that f(R) theories, in general, can ex- 
plain the evolution of the Universe, from some matter 
dominated epoch to the present, late-time self accelerat- 
ing era [12]. One can also generalize f(R) type gravity 
theories by including the function f(R) in the bulk ac- 
tion of the branc- world theories [13|. For this type of 
gen eralized brane world theories it has been shown in 
[14j that they can describe a Universe which starts from 
a matter-dominated era and ends in an accelerated ex- 
panding phase. The classical tests of this theory were 
also considered in [15| . 

A new class of modified theories of gravity was pro- 
posed recently, consisting of the superposition of the met- 
ric Einstein-Hilbcrt Lagrangian with a f(lZ) term con- 
structed a la Palatini [161 ] ■ The dynamically equivalent 
scalar-tensor representation of the theory was also for- 
mulated, and it was shown that even if the scalar field 
is very light, the theory passes the Solar System obser- 
vational constraints. Therefore the theory predicts the 
existence of a long-range scalar field, modifying the cos- 
mological [13, Q3 and galactic dynamics [19| . An explicit 
model that passes the local tests and leads to cosmic ac- 
celeration was also obtained. The cosmological applica- 
tions of a hybrid mctric-Palatini gravitational theory was 
also explored, and cosmological solutions coming from 
the scalar-tensor representation were presented in [l8| . 
Several classes of dynamical cosmological solutions, de- 
pending on the functional form of the effective scalar field 
potential, describing both accelerating and decelerating 
Universes were explicitly found. Furthermore, the cosmo- 
logical perturbation equations were derived and applied 
to uncover the nature of the propagating scalar degree of 



freedom and the signatures these models predict in the 
large-scale structure. Applications to wormhole physics 
were also explored [20I ]. 

Another interesting fact of f(R) gravitational theories 
is that they are equivalent to Brans-Dicke theories with 
a specific u parameter (9j] . This suggests that the theory 
describes the non-minimal coupling between matter and 
geometry in the Einstein frame. It also provides a moti- 
vation to consider non-minimal coupling between matter 
and geometry in a more general manner at the action 
level. In the Einstein-Hilbert action, which has a sim- 
ple additive structure in terms of R and L m , geometry 
and matter appear at two very different conceptual lev- 
els, without any interaction between them. However, the 
idea that the gravitational action may not be additive in 
matter and geometry cannot be rejected a priori. One of 
the first efforts in this direction was made in 21] where, 
based on very general physical arguments, a framework 
was suggested in which terms with non-minimal coupling 
between matter and geometry should be expected in the 
action. As a consequence, a general action for the grav- 
itational field would require a general coupling between 
the Ricci scalar and the matter Lagrangian. 

In this context, a maximal extension of the Einstein- 



Hilbcrt Lagrangian was introduced in 22| , where the La- 
grangian of the gravitational field was considered to be 
a general function of R and L m and therefore this the- 
ory came to be known as the f(R,L m ) gravity theory. 
In theories with non-minimal geometry-matter coupling 
there exist an extra force, which arises from the interac- 
tion between matter and geometry, as initially suggested 
in [23[ and [2J| , respectively. This extra force affects the 
motion of a test point particle, causing it to undergo a 
non-geodesic motion [25( . On the other hand it has been 
suggested that the extra force could be ignored if one 
uses a matter Lagrangian of the form L m = p instead of 
L m = —p [26J. However, in [27| it was shown that when 
the particle number is conserved, the Lagrangian of a 
barotropic perfect fluid is C m = — p[c 2 + J P(p)/p 2 dp], 
where p is the rest mass density and P(p) is the pressure. 
In turn, the / (R, L m ) theory was generalized recently by 
considering a gravitational theory with an action given 
by an arbitrary function of the Ricci scalar, the mat- 
ter Lagrangian density, a scalar field and a kinetic term 
constructed from the gradients of the scalar field, respec- 
tively [281 ] . 

Another difficulty of the non- minimal theories is that, 
in general, the equivalence principle is violated. In fact, 
it has been shown that the observational data of the 
Abcll Cluster A586 exhibits evidence of the interaction 
between dark matter and dark energy, and that this in- 
teraction implies a violation of the Equivalence Principle 
[29(. The mass profile in this particular cluster is ap- 
proximately spherical, and it is a relaxed cluster, since 
it has not undergone any important merging process in 
the last few Gyrs. For the Abell Cluster A586 the ki- 
netic energy px and the gravitational potential energy 
pw can be computed. Then the generalized virial theo- 



rem 2px + pw = £,Pw, where £ is a coupling constant, 
allows to estimate the magnitude of the dark energy- 
dark matter interaction, as well as the degree of violation 
of the equivalence principle that should be detectable in 
large scale cluster surveys [29| ■ 

In the f(R,L m ) type theories |22h29I | it is assumed 
that all the properties of the matter are encoded in the 
matter Lagrangian L m . An alternative view would be 
to consider theories in which matter, described by some 
of its thermodynamic parameters, different from the La- 
grangian, couple directly to geometry. In the standard 
ACDM model the cosmological constant is spatially uni- 
form and time independent, as required by the principle 
of general covariance. Physically, it can be interpreted 
as a relativistic ideal fluid obeying the equation of state 
p + p = 0. Moreover, this cosmological fluid obeys an 
equation of continuity that does not depend on the mat- 
ter energy density. Such a form of dark energy is said to 
be non-interacting [30( . An interaction between ordinary 
matter and dark energy can be introduced in the form 
of a time-dependent cosmological constant. However, to 
preserve the general covariance of the field equations, a 
variable cosmological constant must depend only on rel- 
ativistic invariants. The assumption A = A(i?) leads 
to the f(R) class of modified gravity theories. Another 
choice, in which the cosmological constant is a function 
of the trace of the energy-momentum tensor T, was pro- 
posed in [3(|. One advantage of this choice, as compared 
to /(i?)-type gravity theories, is that we do not need to 
transform conformally from one metric to another, and, 
since we use one and the same metric tensor, the problem 
about which metric is physical does not appear. 

Following the initial work done in [30( , the general non- 
minimal coupling between matter and geometry was con- 
sidered in the framework of a Lagrangian of the form 
f(R, T), consisting of an arbitrary function of the Ricci 
scalar and the trace of the energy- momentum tensor [3l| . 
The gravitational field equations in the metric formalism, 
as well as the equations of motion for test particles, which 
follow from the covariant divergence of the stress-energy 
tensor, were obtained. The equations of motion of test 
particles were also obtained from a variational principle. 
The motion of massive test particles is non-geodesic and 
takes place in the presence of an extra force orthogonal 
to the four- velocity. 

The astrophysical and cosmological implications of the 
f(R, T) grav ity theory have been extensively investigated 
recently (32| ■ A reconstruction of the cosmological mod- 
els in f(R,T) gravity was performed in J33(. The dust 
fluid reproduces the ACDM cosmology, the phantom- 
non-phantom era, and the phantom cosmology. The nu- 
merical simulation for the Hubble parameter shows good 
agreement with the BAO observational data for low red- 
shifts z < 2. The study of the evolution of scalar cosmo- 
logical perturbations was performed [3J], by assuming 
a specific model that guarantee the standard continuity 
equation. The complete set of differential equations for 
the matter density perturbations was obtained and it was 



shown that for general f(R,T) Lagrangians the quasi- 
static approximation leads to very different results as 
compared to the ones derived in the frame of the ACDM 
model. For sub-Hubble modes, the density contrast 
obeys a second order differential equation, with explicit 
wave-number dependence, and subsequent strong diver- 
gences on the cosmological evolution of the perturba- 
tions. A comparison of these results with the usual quasi- 
static approximation in general relativity shows that the 
density contrast quantities evolve very differently. There 
is also a difference in the linear regime between these 
theories. The results obtained in [34[ for f(R,T) gravity 
are in contradiction with the usually assumed behavior of 
the density contrast, and imposes strong limitations on 
the viability of the f(R,T) = f 1 (R) + f 2 (T) type models. 
The growth of the scalar perturbations in the sub- Hubble 
limit, for this model, is scale-dependent. However, one 
should emphasize that the observational data provided 
by the Planck satellite [35[ show clear evidence of the 
scale dependence of the CMB power spectrum. On the 
other hand it seems that Lagrangians of the form f(R, T) 
cannot lead in general to the standard energy-momentum 
conservation equations [3l|, |34j • 

However, the f(R,L m ) or f(R,T) type theories are 
not the most general Lagrangians describing the non- 
minimal coupling between matter and geometry. For ex- 
ample, one may generalize the above modified theories of 
gravity by introducing a term R^JT^" in the Lagrangian. 
Indeed, examples of such cou plin gs can be found in the 
Einstcin-Born-Infcld theories [36[ when one expands the 
square root in the Lagrangian. An interesting difference 
in f(R,T) gravity and in an inclusion of the i? /i „T A "' 
term, is that in considering a traceless energy-momentum 
tensor, i.e., T = 0, the field equations of f(R,T) grav- 
ity reduces to those of f(R) gravity theories. However, 
considering the presence of the R^JT^ coupling term 
still entails a non-minimal coupling to the electromag- 
netic field. 

It is the purpose of this work to consider an exten- 
sion of the f(R, T) gravity theory by also taking into ac- 
count a possible coupling between the energy-momentum 
tensor of ordinary matter, !),„, and the Ricci curvature 
tensor Rp, v . Therefore we propose to describe the grav- 
itational field by means of a Lagrangian of the form 
/ (R,T, R,,vT^ u ) (a similar approach is carried out in 
[371 ]. but in a different setting), where / is an arbitrary 
function in the arguments R, T, and i? AU/ T A " y , respec- 
tively. We obtain the gravitational field equations for 
this theory, and formulate them as an effective Einstein 
field equation. The equation of motion of massive test 
particles is also obtained from the field equations. In 
this type of theories the energy-momentum tensor is gen- 
erally non-conserved. In order to study the Newtonian 
limit of the theory we derive the equation of motion from 
a variational principle. An important requirement for 
any generalized gravity theory, besides passing the Solar 
System tests, is its stability. Thus, we analyze in detail 
the so-called Dolgov-Kawasaki instability, obtaining the 



stability conditions for the theory. An interesting ques- 
tion is the possibility of the conservation of the energy- 
momentum tensor in such theories. We impose the con- 
servation of the energy-momentum tensor by employing 
a Lagrange multiplier method. The gravitational equa- 
tions with energy-momentum conservation are derived 
from an action with the Lagrange multiplier, multiplying 
the energy-momentum tensor, included. The cosmolog- 
ical implications of the theory are investigated for both 
the conservative and non-conservative cases, and several 
classes of analytical and numerical solutions are obtained. 
The present paper is organized as follows. The gravita- 
tional field equations of / (R, T, R^T^) gravity theory 
are derived in Section HH and the equations of motion 
of massive test particles are obtained in Section IIIII The 
Newtonian limit of the theory is studied in Section IIV1 
where in particular, we obtain the generalized Poisson 
equation. In Section [Vl the Dolgov-Kawasaki instability 
in the / (R,T, R^JT^) gravity theory is further investi- 
gated. In Section lVTl the field equations with a conserved 
energy-momentum tensor are obtained via the Lagrange 



multiplier method. In Section IVIII the cosmological im- 
plications of the theory are investigated. Wc discuss and 
conclude our results in Section FVIIII We work in a system 
of units with c = 1 . 



II. THE FIELD EQUATIONS OF THE 

/ (i?, T, i^T"") GRAVITY THEORY 

We consider that the non-minimal coupling between 
matter and geometry can be described by the follow- 
ing action, containing, in addition to the Ricci scalar 
R and the trace of the energy-momentum tensor T, an 
explicit first order coupling between the matter energy- 
momentum T^ v and the Ricci tensor, respectively, 



S = 



16ttG 



^ f ^x^f~gl (R, T, R^T^) + f d A xyf~gL m , 

(1) 

where L m is the Lagrangian density of the matter sector, 
and the matter energy-momentum tensor T^ v is defined 
as 



T — — - 



2 5_W=gL« 
f=g Sgv 



Qfjhu^m & 



. dL m 



(2) 



In the second equality we have assumed that the La- 
grangian is a function of the metric and not its deriva- 
tives. The only requirement imposed on the function 
/ (R, T, R^T^) is that it is an arbitrary analytical func- 
tion in all arguments. 

The variation of the action with respect to the metric 



can be written as 

1 



SS = 



16ttG 



rf 4 



xv -g 



f R SR + f T ST + f RT 5{R^T^ 



* v 9 



r .^v 



(3) 



where we have defined /#, fx and Jb,t as the derivatives 
of / with respect to R, T and R iiV T ilv , respectively. 
The variation of the Ricci scalar is provided by 

5R = R^SfT + g^(V x 5T\ v - V„5T\ X ), (4) 

and variation of the Christoffcl symbol is given by 

<$rV = 2-9 A "( V A^ + V^Sop - V a 8g^). (5) 

Hence, variation of the Ricci scalar can be written as 

5R = R^SgT + g^DSg^ - V^VJg^. (6) 

In addition, variation of the trace of the energy- 
momentum tensor yields the following 



& T _ ST a/3 g/3 ,rp 

,9 T J-fxv 



SgV 5g»» ■ 



(7) 



Using Eq. (J5|), for variation of the energy- momentum ten- 
sor wc find 



ST 



a/3 



= dg a/: 
5g^ dgV" 



L, 



9a/3 



dL Ti 



d 2 L r . 



dg^ v dg^dg a P 



dgo£ 1 t I t 9 d 2 ^m 

daw m 2 Sal39flu 2 9c 'P ^ dq^dq^' 



(8) 



In order to obtain the variation of i? // „T A " / , we note 
that variation of the Ricci tensor is given by 

sr^ = VxSr\, - vjr\ x 

= V Q V ( ^ 5 , )Q - ±n6g M u + g a pV^„6g a ?. (9) 
Therefore wc obtain 



S(R^T^) = T^V A V (Q % )A - -T ap DSg aP 



aT al3 



+T ap g^V a V p 5g^ + R^-^—Sg^. (10) 

Combining all these results, the gravitational field equa- 
tions corresponding to action Eq. (JTjl takes the following 
form 



if Ft, — JRT^m)Gy,v 



n/fl + \r!b, - \f + hUa + \v a v p {.f RT T a P) 



g^u - V^V^/i? + -□ [farT^u) 



2f RT R a(lx T v f - V Q V (M T a v) f RT - ( / T + i/icrfl + 8^g) T m „ - 2 (/ T .^ + / flT i?^" 



dgvdg"? 



0. (11) 



r 



The trace of the gravitational field equation, Eq. fjlljl , 
is obtained as 



and 



3Df R + l -U (f RT T) + V«V^ (f RT T ap ) + Rf R - Tf T 
-\RTf R T + 2R aP T a ?f RT + Rf RT L m + if T L m - 2/ 



-SttGT - 2g^ (grtfr + R afi f RT ) 



d 2 L r 



QgltvQgap 



0.(12) 



T eff _ 1 

^ v fn. - f R rL r 



1 



V^Vvfn - V a f RT V a T^ 



^IrtOT^ — 2f RT R a ^TJ^ + V Q V(^ T^fnT 



2 {f T g afj + f RT R 



a/3\ 



d 2 L r 



g g ^g ga /3 I 



(17) 



The second derivative of the matter Lagrangian with 
respect to the metric is non-zero if the matter Lagrangian 
is the second or of higher order in the metric. Thus, 
for a perfect fluid with L m = —p, or a scalar field with 
L m = —d^(j)d^4 , /2 1 this term can be dropped. However, 
for instance, considering the Maxwell field, we have L m = 
—F^F^/4, and this term results in 



d 2 L r . 



dgvdg a P 



' n? fia-t'uP-, 



(13) 



thus giving a non-zero contribution to the field equations. 
In the framework of / (R, L m ) theories it has been shown 
in [26| that for a matter source in the form of a perfect 
fluid, for a non-minimally coupled Ricci scalar and mat- 
ter Lagrangian in the form L m = p, the extra force van- 
ishes in the case of dust. However, in the present case, we 
will see that even with this choice, the extra force does 
not vanish in general. 

In analogy with the standard Einstein field equation 
one can write the gravitational field equation (|lll) as 

G^ v = 8irG eff T^ - Kff9 M u + T e JJ , (14) 

where we have defined the effective gravitational coupling 
G e ff, the effective cosmological constant A e //, and an 
effective energy-momentum tensor T*l* as 



respectively. In general G e // and h. e ff are not constants, 
and they depend on the specific model considered. 

It is worth mentioning the main differences between 
the present theory to that presented in [3l|. In particu- 
lar, when assuming a traceless energy-momentum tensor, 
T = 0. For instance when the electromagnetic field is in- 
volved, the gravitational field equations for the f(R,T) 
theory reduce to that of the field equations for f(R) grav- 
ity and all non-minimal couplings of gravity to the matter 
field vanish. In contrast, the theory outlined in this work 
still has a non-minimal coupling to the electromagnetic 
field via the R^T^ coupling term in the action, which 
is non-zero in general. It is interesting to write down the 
field equations for this specific form, which are given by 



1 



f R + -J R TF a0 F^\G^ 



-\f+ l -V a Ve{f RT T^) 



a.f R + -Rf R 



g^y - V M V„/j 



1 



□ UhtT^v) + 2f RT R a ( tl T v " - V Q V( M T a v ^f RT 



\fiwR + 8irG ) T^ + f RT R ati F m F v& = 0. (18) 



G + ±(f T + \fwrR - |D/gr) , . 

G eff = 7 7 — j > (15) 

J R — JUTJ-'m 



III. EQUATION OF MOTION OF THE 

MASSIVE TEST PARTICLES IN THE 

f {R,T,R liv T>"') GRAVITY THEORY 



A 



eff 



2Df R + Rf R -f + 2f T L m + V a V (f RT T a P) 



2(f R — f R rL r 



(16) 



The covariant divergence of the energy-momentum 
tensor can be obtained by taking the divergence of the 
gravitational field equation, Eq. (fTTj) . which takes the 



following form 



V' J T = 

"" (l + Rfa R + 2fa) 



V M (f RT R afi T au ) 



V„ (L rn f T ) - - [ f RT R pa + fag pa ) vjr 9 " 



~ G^V (f RT L m ) - - [W (Rf RT ) + 2V7t] V y 

(19) 

where we have assumed that d 2 L m jdg^ v dg a ^ = 0, and 
we have used the mathematical identities 

V^ (f R R^ + Df R9fW - ±f giu , - V„V„/.r 



fTVvT + f K rVu(Rp*T'") 



(20) 



^^/jr;«5[;p;cr] — ± fiT.ati g pa + J^aT;S^ ppa > -LfJ.a:5li- t pen 

(21) 
and 

[n,V u }T = R^WT, (22) 

respectively. 

In order to find the equation of motion for a massive 
test particle we start with the energy-momentum tensor 
of the perfect fluid, given by 

Tfj, v = pg pv + {p + p)u p u u , (23) 

where u M is the four-velocity of the particle. Taking the 
divergence of the Eq. (|23)l , and by introducing the pro- 
jection operator h^, defined as 



Q\lv \ ^p^v , 



(24) 



we obtain 

+ (p+p)(«"V/ + u"V/), (25) 

Multiplying the above equation with h x one finds 

h * VfiT n» = (j> + p )^v M u A + h vX V v p, (26) 

where we have used the identity u^V^u^ = 0. The equa- 
tion of motion for a massive test particle with the matter 
Lagrangian L m = p, then takes the form 



d 2 x x 
ds 2 



T\y i u» = f\ 



(27) 



where we have used the definition of the covariant deriva- 
tive to obtain the LHS of the above equation. The extra 
force acting on the test particles is given by 



f 



1 



(fa + RfRT) V„p - (1 + 3/ r ) V„p 



P + P 

(p + p)f 7 p h av ) 

- f R TR* P h ap V u { P + p) 



1 + 2 fa + Rf RT 



(28) 



Contrary to the nonminimal coupling presented in [23J, 
and as can be seen from the above equations, the extra 
force does not vanish even with the Lagrangian L m = p. 
The extra-force is perpendicular to the four-velocity, 
satisfying the relation f x u\ = 0. In the absence of 
any coupling between matter and geometry, with fa = 
fur = 0, the extra- force takes the usual form of the 
standard general rclativistic fluid motion, i.e., / A = 
-h Xv V u p/ (p + p). In the case of / (R, T, R^T^) grav- 
ity theories, there is an explicit dependence of the extra- 
force on the Ricci tensor R ap , which makes the devia- 
tion from the geodesic motion more important for regions 
with strong gravitational fields. 



IV. THE NEWTONIAN LIMIT OF 
f{R,T,R pl/ T^) GRAVITY 

Lets us now consider the Newtonian limit of the theory. 
Using the weak field and slow motion approximation, we 
derive the equation of motion of massive test particles 
in a weak gravitational field as well as the generalized 
Poisson equation satisfied by the Newtonian potential <\>. 



A. The equation of motion of massive test particles 

In order to obtain the Newtonian limit, we show first 
that the equation of motion, Eg. (|2T1). can be derived 
from a variational principle [2J, |31| . To this end, we 
assume that one can represent the extra force formally 

as 



f X = {g vX + U v U X )^ v \n^Q, 



(29) 



where Q is a function to be determined from the vari- 
ational principle. With this assumption, one can prove 
that the equation of motion Eq. ([27)) can be obtained by 
varying the action [22( 



L p ds 



Q\Jg^u^u v ds, 



(30) 



where S p and L p are the action and Lagrangian density 
of the test particle respectively, provided that Q is not 
an explicit function of u^ . When \[Q — > 1 , we obtain the 
variational principle for the standard general relativistic 
motion for a massive test particle. 

In order to obtain the function Q for Eq. (|27[) in the 
Newtonian limit, we assume that the density of the physi- 
cal system is small and therefore the pressure satisfies the 
condition p<p. Hence the energy-momentum tensor of 
the system can be taken as the energy-momentum tensor 
of pressureless dust. Moreover, by considering the limit- 
ing case of small velocities, we can take the four- velocity 



in the form u^ = Sq / y 
tives of h pv in Eq. 
form 



<?oo and drop the covariant deriva- 
Thcrcforc Eq. (gTJ takes the 



f X = -h x »\7„p, 
P 



(31) 



where 



F = 



fT + fR.T(R-R a ph a P) 
1 + 2/ T + RIbt 



(32) 



In the Newtonian limit, one can expand the quantity F 
as a function of the energy density as 



F(p) = F(p )- 



dF 
dp 



(p-po)=F +F 1 (p- Po ). (33) 



Substituting F to the equation (|31|) and comparing the 
result with equation (f2U)) one can read off \[Q for small 
p as 



(34) 



yfQ-l + Fyp- 



where a = Fq — Fip^. We have therefore obtained y<9 
in the case of dust as an explicit function of the energy 
density p. We may now proceed to study the Newto- 
nian limit of the theory by using the variational principle 
Eq. ([30]), and also Eq. fl34|). In the weak field limit the 
interval ds for dust moving in a gravitational field is 



ds ps y/l + 24>-v 2 dt ps 1 



- - 1 *, 



(35) 



where is the Newtonian potential and u is the three- 
dimensional velocity of the fluid. The equation of motion 
of the fluid to first order approximation can be obtained 
from the variational principle 



i + u( P ) + 4>-- 



dt = Q. 



The total acceleration of the system, a, is given as 
a = -V(/> - VU(p) = a N + Be, 



(36) 



(37) 



where cln = — V0 is the Newtonian acceleration, and the 
supplementary acceleration, induced by the geometry- 
matter coupling, is 

a E (p) = -VU{p) = [F x + (F - Fip )p Q_1 ] Vp. (38) 



The acceleration given by Eq. (1351) is due to the modifi- 
cation of the gravitational action. In our case, there is no 
hydrodynamical acceleration a p term in the total accel- 
eration, because of our assumption that the fluid is pres- 
sureless. However, such an acceleration does exist in the 
general case. We see from Eq. (f3"5|) that the extra accel- 
eration cle is essentially due to the non-minimal coupling 
between matter and geometry. The extra-acceleration is 
proportional to the gradient of the energy density of the 
fluid. Therefore, for a constant energy density source and 
a prcssureless fluid, the extra acceleration vanishes. 



one has R = — 2i?oo = — 2V 2 0, where is the Newtonian 
potential which appears in the (00) component of the 
metric goo = — (1 + 20), one can compute the individual 
terms in the trace equation (fT2"j) as 



na/3 



R a pT ap ~ pW 



□/r ~ V 2 /i? + V/r • V0, 



(39) 



(40) 



and 



V^VM/htT"' 3 ) ~ V(pf RT ) ■ V0 + pf RT V'4>, (41) 

respectively. 

Substituting the above expressions into Eq. (fT2"|) and 
rearranging terms, we obtain the generalized Poisson 
equation as 



1 



2(/r - 2pf RT ) 

2f + V(3f R + pf RT )-V<f> 



8irGp + 3W 2 f R - 3pf T 



(42) 



As can be seen, the generalized Poisson equation is mod- 
ified by the addition of gradient of the <j) field to the 
equation. 



V. THE DOLGOV-KAWASAKI INSTABILITY 

IN f(R,T,R^T^) GRAVITY 

Beside consistency with the Solar System tests, any 
gravitational theory should be stable against classical and 
quantum fluctuations. One of the important instabilities 
of modified theories of gravity is the Dolgov-Kawasaki 
instability [38l . |39| , which we shall discuss in the present 
Section. 

Let us assume that, in order to be consistent with the 
Solar System tests, the Lagrangian can be written as 

f(R, T, R^T^) = R + e$ (R, T, R^T"") , (43) 

where e is a small parameter. Following [38(, we expand 
the space-time quantities around a constant curvature 
background with geometrical and physical parameters 
{v^,Ro, T° v , T , L ) , so that 



R/J.U = -rRoViu' + Rfiv> R — Ro + Ri, (44) 



T^^Tl + Tl, T = T + T 1 , (45) 



B. The generalized Poisson equation 

To obtain the Poisson equation we assume that the 
matter content of the self-gravitating system is repre- 
sented by dust. Also, noting that in the Newtonian limit 



L m = L + L l7 (46) 

where we have locally expanded the metric tensor as 

SV = Vt^y + hfiv (47) 



We note that in the above equations we have really two 
types of approximations, as mentioned in [39(. The first 
is an adiabatic expansion around a constant curvature 
space, which is justified on the time-scales much shorter 
than the Hubble time. The second approximation is a 
local expansion in the small regions of space-time, which 
are locally flat. These approximations have been used ex- 
tensively in f(R) gravity theories, [3a,|39|. The function 
f{R,T,R liV T' lv ) can be expanded as 



f(R,T,R liU T^) = R +R 1 + e 



$(0) + $ R (0)Ri 



+ a> T (0)T! + $ RT (0) f -RoT 1 + Rljrp 

= Rq + e$(0) + [l + 6*^(0)] iJi + # (1) , (48) 



where (0) denotes the computation of the function at the 
background level, and for simplicity we have defined the 
first order quantity H^ as 



H^=e 



$r(0)T! + $ RT (0) -RoT 1 + R 1 ^ 



1 



nllV 



We then obtain 



f R = 1 + e$ fl (0) + e$*,fl(0)fli + H R > 



(i) 



(i) 



f T = e$r(0) + e$ T ,fl(0)fli + H% 
far = ^ RT (0) + e$ flT ,R(0)iii + H { R %. 



(49) 

(50) 
(51) 
(52) 



The trace equation (TT21 can then be expanded to first 
order to obtain 



3e$ fl ,a(0) + ieT $ Rr ,R(0) ) ORi + eT o Q/3 $ flT , fl (0)V Q V /J i?i 



/ fl (0) + ei? o $fl,i?(0) - eT o $ Tji? (0) - -T / flT (0) 



-^eRoT o <S> RT , R (0) + ^eR o T o <P RT , R (0) + eRoL $ RTiR (0) + f RT {0)L o + 4ei $T,i?(0) - 2 - e$ H (0) 

+3Di/j ? 1) + i/^OpTi + ir D^ + /i ?T (0)V Q V /3 T 1 a ' 3 + T^V^pH^l + RoH R 1] - Ti/ T (0) 

-TqF^ + 2fl* 1/ T '"7 icr (0) + RofMO)^ + RoLoH™ + 4/ r (0)ii + 4L ff£ ) 

-2# « + SttGT! - VY* [/ r (0) + ^o/ht(0)] ^g^ = 0. 

In the limit considered, one may write □ = —d 2 + V 2 , thus obtaining 

T^V a VpRi = T^R, + TtfdidjRt. 
One can then rewrite the above equation as 



Ri 



where wc have defined 



V eff 



R X + V^jfViVjR! + mlfjRi = H eff , 



(3e$fl, fl (0) + \eT $> RT , R (0))5v + eT^ RT , R (0) 
T 00 - 3e$ R)fl (0) - £eT **r lfl (0) 



and 



H eff = [3e$ flifl (0) + ^T o $ fl T,fl(0) - T 00 ] '{sOff^ + ^(0)0^ + \t UH r % + f RT (O)V a V Tf 
+T O!/3 V Q V /3 #W + ifcffW - ri/ r (0) - T H£ } + 2R 1 T^fsriO) + Rof«r{0)Li + R L H R ^ 



-4/ T (0)ii + 4L H^ ] - 2H {1) + 87rGTi - 2rf v r\ 



a/3 



M0) + ^RofRT(0) 



d 2 L x 



QgllvQgOlp 



respectively, and we have introduced the effective mass as 
(T°° 



m eff 



-^T )f RT , R (0)-3f RR (0) 



eR ® R , R (ty - eT $ T ,fl(0) - -eR o T o <S> RT , R (0) 



-\eT^ RT {Q) + ^ei? To$fl T ,fl(0) + eR o L o $ RT . R {0) + e<P RT {0)L o + 4eL $ T ,fi(0) - 1 - e$ fl (0) 



(53) 

(54) 
(55) 

(56) 



(57) 



(58) 



r 



The dominant term in the above expression is 1/[3/hr(0) + {\Tq — TQ°)f R T. R (0)] , an d therefore the 



condition to avoid the Dolgov-Kawasaki instability is 

3/flfl(0) - (p - ^Toj /« T ,fl(0) > 0, (59) 

where po is the background energy density of the matter 
[381 ] . We note that due to the above expression, condition 
for the stability does not depend on the derivative of 
the function / with respect to T. So, the DK stability 
condition for the case of f(R,T) gravity is the same as 
f(R) gravity. However, the condition is modified in the 
case of f(R,T,R^T^). 



a barotropic fluid, the function J2{T) can be uniquely 
determined as /2(T) ~ T 1 / 2 . In the following we investi- 
gate the energy conservation in / (R, T, R IJ , I/ T flv ) gravity. 
In order to impose the matter energy-momentum ten- 
sor conservation, one can use the Lagrange multiplier 
method J4Jj. To do effect, let us consider the modified 
action 



S = 



1 



j4 i 

a Xy/—g 



16ttG 



/(R^R^T^ + X^T^ 

(60) 



VI. f{R,T,R^T^) GRAVITY THEORIES WITH 
ENERGY-MOMENTUM CONSERVATION 

The general non-minimal coupling between matter and 
geometry leads to the important consequence that the 
matter energy-momentum tensor is not conserved. In 
Section irm we have shown that this property of the grav- 
itational theory determines the appearance of the extra 
force. However, the energy non-conservation can be in- 
terpreted as a shortcoming of these types of theories. In 
the framework of the f(R,T) theory, models with en- 
ergy conservation have been investigated in [34j . By as- 
suming a specific additive form for the function f(R,T), 
f(R, T) = f\{R) + f2{T), and by imposing the condition 
of the energy conservation, and under the assumption of 



where we have introduced the vector Lagrange multiplier 
A M . The variations of the first and the third terms are 
similar to those computed in Section ITT] The variation of 
the second term with respect to the metric is given by 



{(^"VT^) 



A a (V M T ar , - -V^T aP g^)dg^ + X^VST^ 



(61) 



where the variation of the energy-momentum tensor is 
given in Eq. (|8]). Combining the above results with the 
calculations of Section UH we obtain the field equations 
together with the energy-momentum conservation as 



{Jr. - JKrRm) G^ + 



af R + \rJr - \.f + f T L m + iv Q v^ {f RT T^) 

1 



SV - V m V„/h + -D(/ flT T /Jiy ) 



8 2 L r 



+2f RT R a ( p T l/ f - V Q V (/i [T a v) f RT \ - (It + ^IrtR + 8ttG)T^ - 2(f T g a ^ + f RT R ap ) Q ^ v ™ afs 



d 2 L r 



\\ a V p T ap9tlv + A P V (M T; - V (fl \ v) L m - \v a \ a (L m g^ - T„ v ) + 2V {a \ 0) ^- 



Now, variation with respect to the vector X p results in 

V%„ = 0, 



dg^ w dg a P 



0. 



(62) 



(63) 



which is the conservation of the energy-momentum tensor. Therefore Eqs. (|62p and (|63p provide the basic equations 
of the / (R,T, R^T^") gravity theory with energy conservation. The gravitational field equations explicitly depend 
on the Lagrange multiplier A M . The field equations in the case of the matter Lagrangian L m = —p or L m = p which 
leads to d 2 L m /dg pu dg a ^ = take a simpler form. For L m = —p we obtain 

Ur + p!rt)G^ + pf R + l -R.f R . - \.f - pfr + \T aP V a V p f RT \g^ - V M V„/ fl + l -U(.f RT T^) + 2f RT R Q ^T„f 

- V Q V (M [T a v) f RT \ - (/t + \f RT R + 87rG)T^ + X P V (li T p u) + pV^X v) + \{pg^ v + T^)V Q A« = 0, 

(64) 



r 



where the conservation of the energy-momentum tensor is taken into account. 
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In the case of the electromagnetic field, because the 
trace of the energy-momentum vanishes, we have /t = 0, 



and using Eq. (fT5|) we find the field equations 



(!r + \f 2 !et)G^ + [□/* + \Rf R - i/ - pf T + \T^V a V f RT ]9^ - V M V„/ fl + ]p{f RT T^) + 2f RT R a{ ^ 



V Q V (Al [T a v) f RT ] - [\f RT R + 8ttG)T^ + f RT R ap F^F H 



WfrK: 



-F V^Xu) 



+ \F 2 ^ a \ a - V (a A«F^i^ = 0, 



(65) 



r 



iF al3 and F^ 



F F a 



where we have defined F = F a p 

and used the conservation of the energy-momentum ten 

sor. 



VII. COSMOLOGICAL APPLICATIONS OF 

/ (R, T, R^T"") GRAVITY 

Let us now consider some examples of cosmological 
solutions of the theory. In order to obtain explicit re- 
sults and as a first step, one has to fix the functional 
form of the function / (R,T,R^ V T^ U ). In the follow- 
ing we consider three specific choices for /, namely 
/ = R + aR„ v T^, f = R + aR^T^ + (3VT and 
/ = R+aRR^T^, where a, (3 = constant, respectively 
We analyze the evolution and dynamics of the Universe 
for the above with and without energy conservation. In 
all cases we assume that the Universe is isotropic and 
homogeneous, with the matter content described by the 
energy density p, and thermodynamic pressure p with 
the matter Lagrangian as L m = —p. The geometry of 
the space-time is described by the Fricdmann-Lemaitre- 
Robertson- Walker (FLRW) metric, given by 



ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), 



(66) 



where a(t) is the scale factor of the Universe. We define 
the Hubble parameter as H = a/a, and we describe the 
accelerated expansion of the Universe through the values 
of the deceleration parameter q, defined as 



Q 



d 1 
dtH 



1. 



(67) 



If q < 0, the expansion of the Universe is accelerating, 
while positive values of q, q > 0, describe decelerating 
evolutions. 



A. Specific case I: f = R + aRy, v T» v 

Let us first consider the simplest case, in which the 
interaction between matter and geometry takes place 
only via the coupling between the energy-momentum 
and Ricci tensors. This simple case can also show the 



main differences of the present theory with the so-called 
f(R, T) gravity theory [3lj . The gravitational field equa- 
tions for this form of / are given by 



G 



{lis 



^^(T(fl^ t/) Q-^p<T-^ Qj-LV c.tx^-^lV 



- (2V ff V ( „T; ) - DT^ - V Q V /3 T Q 'V 



G iwi^r. 



2R 



a/3 



d 2 L r 



dg^ v dg a P 



8irGT^ = 0. (68) 



The effective gravitational coupling, the effective cos- 
mological constant, and the effective energy-momentum 
tensor arc given for this choice of / by 



Geff 
Ae// 



16ttG + aR 
167r(l - aL m y 

97i °Lr d VcVp-Rcp)^, 
2(1 - aL m ) 



(69) 
(70) 



and 



T eff = a 

» v " 2(1 -aL m ) 



5/X/3 V a V^ + gp^aVfj, 



-g^agyp^ - 4i? Q(M g 



„)n 



T 



a/3 



(71) 



For the case of the FLRW metric the independent cos- 
mological field equations are 

3H 2 = ^—s+l 1 -?L-H(p-p), (72) 



1 — ap 2 1 — ap 



and 



2H + 3H 2 



2o 



1 + ap 



Hp 



tip 



1 a 



1 + ap 2 1 + ap 



(P-P) 



(73) 

respectively, where we have denoted n = &t:G for sim- 
plicity. When a = we recover the standard Fricdmann 
equations. To remove the under determinacy of the field 
equations, we must impose an equation of state for the 
cosmological matter, p = p(p). A standard form of the 
cosmological matter equation of state is p = up, where 
lo = constant, and < lo < 1. 
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High cosmological density limit of the field equations 



C\ must satisfy the constraint 



We shall first consider the high energy density limit of 
the system of modified cosmological equations (|?2")) and 
([73)) . Moreover, we assume that the constant a is small, 
so that ap <C 1, and ap <C 1, respectively. In the high- 
energy limit, p = p, and Eqs. (J73J) and (|73|) take the 
approximate form 



3H 2 
2H + 3H 2 



K,p, 

—up + 2aH p. 



(74) 
(75) 



The time evolution of the Hubble parameter is de- 
scribed by the equation 



1-55! | 

K 



3H 2 = 0, 



(76) 



and hence for this model the evolution of the Hubble 
parameter is given by 



A /(Ci+3 K t)2-24a« + Ci+3«t 

#W = 7^ > (<<) 

via 

where C\ > is an arbitrary integration constant. The 
scale factor of the Universe is given by 



exp 



a{t) = C 2 



(Ci+3ret)^/(Ci+3Kt) 2 -24aK+9Kt 2 +6KCit 
72ok 



^V(Ci + 3Ki )' 



24«k + Ci + 3k£ 



(78) 

where Ci > is an arbitrary integration constant. In 
order to have a physical solution the integration constant 



d > V24k«. 



(79) 



The values of the integration constant can be deter- 
mined from the condition H(0) = Hq, and a(0) = ao, 
where i?o & n d ao are the initial values of the Hubble pa- 
rameter and of the scale factor of the Universe, respec- 
tively. This condition immediately provides for C\ the 
following value 



Ci 



6aiJg + « 



(80) 



For the integration constant Ci we obtain 



C 2 = 



ao 



1 



'(/c-6affg)' 



6aHn 



Hn 



x exp 



72ai? K 



(81) 



In the small time limit, the scale factor can be repre- 
sented by 



a{t) w oo 1 + 1 r^—t 
bn a 



The deceleration parameter is obtained as 



(82) 



J 



q{t) 







36a Hqk 




(6aH->+3H Kt+K) 2 

V HS 


- 24a« 


&aH* + H ^ 6aH2 ° +3 H H r t+K T- 


- 24a« + 3H nt + k 



(83) 



and it can be represented in a form of a power series as 



2. The case of dust matter 



?(*) 



WaH$ 



6H k 2 



6aH 2 { K _ Q aH 2 ) 3 



Next we consider the case of low density cosmological 

matter, with p = 0. Moreover, we assume again that 

(84) the condition ap<l holds. Then the gravitational field 

equations, Eqs. (f72"j) and (J73J), corresponding to a FLRW 

Universe, take the approximate form 



For small values of time, if 2AaH^ Ck^w-I, and the 
Universe starts its expansion from a de Sitter like phase, 
entering, after a finite time interval, into a decelerating 
phase. On the other hand, if n > QolHq, q < — 1, and 
the non-singular Universe experiences an initial super- 
accelerating phase. 



3H 2 = np+ -aHp, 
2H + 3H 2 = 2aHp+-ap. 



(85) 



(86) 



First we consider the matter dominated phase of the 
model, in which the non-accelerating expansion of the 
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Universe can be described by a power law form of the 
scale factor, so that a = t m , m = constant, and H = 
m/t, respectively. The deceleration parameter is given 
by q = 1/m — 1 Therefore Eq. (|55|) gives for the time 
evolution of the density the equation 



Sam 



p + up — 3- 



2t r ' ' r ' t 2 
with the general solution given by 



(87) 



P(t) 



3poae s» 



Ei 



'(ft) 



Ei 



.'!(> 



3o 



where Ei(z) = — J_ e t dt/t is the exponential integral 
function, and we have used the initial condition p (to) = 
pa. By substituting the expressions of the density and 
of the Hubble parameter into Eq. ([55)1 . to first order, we 
obtain the following constraint on to, 



9to 2 - 10m + 1 
it 2 



0(t 2 



0, 



(89) 



which is (approximately) satisfied if to is given by the 
algebraic equation 9to 2 — IOto +1 = 0, having the solu- 
tions mi = 1, and mi = 1/9, respectively. The decel- 
eration parameters corresponding to these solutions are 
qi = 0, and qi = 8, respectively. Since a value of the 
deceleration parameter of the order of q = 8 seems to 
be ruled out by the observations, the physical solution 
has a scale factor a = t, and q = 0. The cosmological 
solutions with zero value of the deceleration parameter 
are called marginally accelerating, and they describe the 
pre-accelerating phase of the cosmic expansion. 

Now we look for a dc Sitter-type solution of the field 
equations for the pressureless matter, Eqs. ([85]) and (|86|) . 
by taking H = Hq = constant. Then it follows that, 
in order to have an accelerated expansion, the matter 
density must satisfy the equation 



2k 
p - H a p -\ p = 0, 



with the general solution given by 



(90) 



p(t) 



e ^H (t-t ) 



fy/apfxa-Hopo) 



\ ^aH 2 - 8, 



yfaH 2 — 8k 



2V& 



(t-to) 



cosh 



y^aHy —8k 



(t-h) 



(91) 



where we have used the initial conditions p(to) = Po, and 
p(to) = Poij respectively. Therefore, in the presence of 
a non-trivial geometry-matter coupling, once the evolu- 
tion of the matter density is given by Eq. (|9Tj) , the time 
evolution of the Universe is of the de Sitter type. 



B. Specific case II: f = R + P^/\T\ + aR^T"" 

In this section, we generalize the previous action by 
adding a term /3vT, P = constant. Such a model, sat- 
isfying the energy conservation, was considered, in the 
framework of the f(R,T) theory, in [3J], where a toy 
model with action given by f(R,T) = R + (3T 1 / 2 was 
investigated. The field equations of the / (R, T, R^T^) 
gravity in this case are 



Gau + ct 



2Ra(tjT v) yRp<?T g^u^—RT^u 



^ (2V CT V ( ,T; ) - DT^ - V Q V /3 T Q 'V 



^fiu-^m 



2R 



a/3 



d 2 L r , 



— 8-kGTuu 



2-y/r 



dg^dg a P 
(L m - T)g^ v - T^v - 2g afi 



/"- 
d 2 L r , 



dg^dg a P 



= 0. 

(92) 



The cosmological equations of this model with the perfect 
fluid matter in the FRW space-time can be written as 

1 



3(1 - ap)H 2 = k P + -aH (p - p) - -£^3^1, (93) 

and 

(1 + ap)(2H + 3H 2 ) = 2aHp - k P 



respectively. 



1 . 2(3p 
+ —a (p — p) , 

2 P) ^H 



High density regime 



(94) 



In the high density cosmological regime the matter 
equation of state is given by the Zeldovich stiff causal 
equation of state, with p — p. Then the field equations 
take the form 



3(1- ap) H 2 = k P - 



(95) 



and 



(1 + ap) (2H + 3H 2 ) = 2aHp -np- V2/3y/p, (96) 



respectively. For a small coupling a, and by assuming 
a/)<l, the field equations reduce to 

3H 2 = K P -V2(3^p-, (97) 

H = aHp-Kp, (98) 

thus giving the evolution equation for the density as 



4 K p 3 / 2 (t) j3Kp(t) - 3V2/3JpTt) 

p(t) = — ¥ == . (99) 

V2/3 {1 - 4ap(t)} + 2KyfpJt) [2ap(t) - 1} 



By neglecting the term ap compared to 1, and by series 
expanding the right hand side of Eq. (IM1) , to first order 
we obtain 



-2V3^o 3/2 1 - 



13 



13 



\/2KyfpJ 

with the general solution given by 



2V2k^pJ ' 

(100) 



/^ 2 



32k 2 



3 tanh 



8 V 3 V S/ 3 






-i 



(101) 



where we have used the initial condition p(0) = po- After 
substituting the density given by Eq. f|101[) into Eq. (p?T|) , 
and performing a series expansion with respect to the 
time, to first order we obtain for the Hubble parameter 



H(t) 



y/2/3 3 - 6V2/3k 2 Po + 8k 3 pI /2 



8k 3/ VV^( k V^-^) 



t. (102) 



In order to obtain a physical solution the parameters 
of the model must satisfy the constraint KyfpE > v2/3. 
For the scale factor of the Universe we obtain 



f J K Po - V2f3y/p~o 
a(t) = a cxp < — 1 

y/2/3 3 - 6V2/3k 2 Po + 8n 3 p 3 Q /2 t 2 



8n 3 / 2 ^J^( K ^-V2(3) 



(103) 



In the high density regime, and in the considered or- 
der of approximation, the expansion of the Universe is 
super-exponential, with the scale factor proportional to 
the exponential of t 2 . 



13 



constant. Then Eq. (|105[) can be immediately integrated, 
to give 



P(t) 



-4Jt (t-t ) JL 

3a 4H 



Poi 



3Hq , 
+ -H J1 (*-*o) 



2 a 



-Po 



Poi 
4ff n 



3 

8a' 



(106) 



where po = p (to), and poi = P (to)- in the limit of large 
time the matter density is linearly increasing in time, 
and hence this model does not have a physical late time 
de Sitter phase. Other types of solutions, including the 
matter dominated phase, can be obtained through the 
detailed numerical study of the system of Eqs. (|104[) and 
(|105p . which will not be performed here. 



C. Specific case III: f = R(l + aR^T^) 

As a third example of a cosmological model we con- 
sider the case in which the function / is given by / = 
R(l + aR^ u T^). The field equations in this case are 
given by 



1 + a{R aP T ap - RL m ) 



G 



tiv 



U(R aP T a P) 



V a V p {RT 



a/3-y 



g^u - aV M V l/ (i? a/ gT 



ad\ 



-aD(RT^) + 2aRR a{p T a y) - aV Q V (M [RT a u) ] 



: aR 2 + 8nG]T, u -2aRR^^^ 



(107) 



In the following, we investigate only the zero pressure 
matter filled Universe. The cosmological gravitational 
field equations are given by 

- 3H 2 + K p + a (l8HHp + \8HHp 

+ 54H 2 Hp - m 2 p + 21H 3 p + 27H i p) = 0, (108) 



and 



2. The pressureless matter fluid case 

In the case of dust, having p = 0, the gravitational 
field equations take the form 

3(1- ap)H 2 = K p+^aHp-^, (104) 
2H + 3H 2 = 2aHp+^ap, (105) 

respectively. We consider the late time expansionary 
phase of the Universe, by assuming a de Sitter type 
form for the scale factor, a(t) — exp(Hot), with H$ = 



-2H - 3H 2 + a 6Hp + l2Hp + 36HHp 
+6Hp + 54HHp + A8H 2 Hp + l5H 2 p + 9H 2 p 
+30H 3 p-9H i p\ =0, (109) 

respectively. The terms proportional to a in the gen- 
eralized Friedmann equations (|108[) and (|109[) play the 
role of an effective supplementary density and pressure, 
which may be responsible for the late time acceleration 
of the Universe. In the following, we look for a de Sit- 
ter type solution of Eqs. (|108|) and (|109[) . assuming that 
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H = Hq = constant. Then the held equations take the 
form 

- 'AH 2 + np + 27H$a (p + H p) = 0, (110) 
-3iJ 2 + H 2 a (9p + 30H o p - 9H 2 p) = 0, (111) 

respectively, leading to the following differential consis- 
tency condition for the matter density p, 



9aHlp + ZaHlp - (36a#o + «) P = 0. 



(112) 



The general solution of Eq. (|112j) is given by 

Pit) 



-iffo(t-to) x \ \/aH Q {H Q p + §p 0l ) 
\ y/U5aH$ + 4k 



x sinh 



y/l45aH$ + 4k 



6y/aH 



(t - t ) 



-po cosh 



6y/aH 



(t - t Q ) 



(113) 



where we have used the initial conditions p (to) = po and 
p(to) = Poii respectively In order that the ordinary 
matter density decays exponentially for t > to, all the 



exponential terms must be negative, which imposes on 
a the constraint a < —k/36Hq. The high energy den- 
sity regime of this model, corresponding to p = p, has 
similar properties with the p = case, that is, it admits 
a de Sitter phase, which can be obtained by analytical 
methods. 



D. / (R, T, RnvT 1 * 1 ') gravity cosmological models 
with conserved energy-momentum tensor 

We now consider cosmological models with a conserved 
energy-momentum tensor. For this case the relevant field 
equations are obtained in Section I VII by using the La- 
grange multiplier method, and arc given by Eqs. (|62[) 
and (|63p . respectively For the isotropic and homoge- 
neous FLRW Universe the energy conservation equation 
becomes 



p + 3H(p + p) =0. 



(114) 



Assuming a barotropic equation of state for the matter of 
the form p(t) = wp(t), lo = constant, and the ansatz A M = 
A(£)<5q f° r the Lagrange multiplier, the gravitational field 
equations with energy conservation are given by 



J 



fn + ^Pfnr) H + 3 ( f R + ^ ±j f" = 2 - pf RT ) H 2 



and 



(w + l)A —Jrt 



, , 3w 2 -w-6 U |9 /, 3w 3 + 3w 2 + 2 \ 2 

/fl + 2 ^ RT J # + 3 ( /« 2 P/j?T J 

jfi - o( w - ^P/ht + r(w + l)pA - (w + 1)p/t - Kujp - -/ = 0. 



p-/fl^ + pA--/ + 87rGp = 0, 

(115) 



(3c 2 -l)/ flT + ^A 



P-Vr}H 

(116) 



r 



respectively, where we have eliminated p from the above 
equations by using the conservation equation Eq. (|114[) . 
As an example for cosmological applications we con- 
sider the case where the function / is given by 



f = R+ CtRavT^, 



(117) 



where a = constant. In this case Eqs. (|115[) and (|116[) 
become 



and 



(3cj 2 -l)ap- 1 



H 2 + 3{uj + l)\pH + pX + k P = 0, 

(118) 



2--(3w 2 -4w-3)p 



H + -auj{uj + ifpH 2 



+ -(5w + 3)X P H + -(1 - w)pA + k(w + l)p = 0, (119) 



respectively 



1. The high energy density phase 

In the high energy density limit we assume that the 
equation of state of the cosmological matter is the stiff 
causal equation of state, with p = p. Then the energy 
conservation equation gives 



Pa 



(120) 



The field equations for the high density phase of the 
evolution of the Universe arc given by 

3(ap-l)H 2 + 6\Hp + pX + Kp = 0, (121) 

and 

(l + ap)H + 9apH 2 + 2\pH + Kp = 0, (122) 
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respectively. By assuming that ap ^> 1, Eqs. (|121[) and 
([T22j) become 



3aH 2 + 6XH + X + « = 0, 
a# + 9aF 2 + 2AF + « = 0. 



(123) 
(124) 



For H = Hq = constant, and A = Ao = constant, and for 
a < 0, Eqs. ([T2"3>(fT2"4]) have the solution 



1 



Ha = l\fS' Xo = i^^- (125) 



Therefore in the / (R, T, R fil/ T til/ ) gravity with energy 
conservation a de Sitter type phase does exist during the 
high density regime of the cosmological evolution of the 
Universe. From the above equations we obtain the re- 
lation between the Lagrange multiplier and the Hubble 
parameter as 



A - 1 K 



(126) 



where we have used the initial condition H (to) = Hq . 
For a < we obtain 



H(t) 



tanh 



tanh 



-Ha 



(t-to) 



(133) 

For a > 0, the time evolution of the scale factor is 
given by 



a(t) = ao cos 2 



^(*-*o)-tan- Uf K H 



(134) 

where ao is an arbitrary constant of integration. For 
a < we have 



a(t) = ao cosh 



(t — to) + tanh 



-Ho 



(135) 

In the opposite limit of small densities ap <C 1, 
Eq. (fnO)) takes the form 



AH + 3H 2 



np 



= 0, 



(136) 



2. The pressureless matter case 

In the case of dust matter, i.e., to = 0, from the con- 
servation of the energy-momentum tensor we obtain the 
density of the Universe as 



Po 



(127) 



The gravitational field equations with the conservation 
of energy-momentum and dust matter take the form 



and 



3(l + -p)H 2 + 3\pH + p\ + Kp = 0, (128) 



2 + ^p\H+^XpH+±pX + K p = Q 1 (129) 



respectively From Eqs. (|128[) and (J129I) we immediately 
obtain 



In the limit of large densities ap ^> 1, Eq. (|130l) be- 
comes 



3aH + 3-H 2 + k = 0, 



with the general solution given for a > by 



(131) 



2k 
Hit) = \ — tan 
V 3a 



tan 






ba 



(132) 



with the general solution given by 
{(3ogflg - K p )f 3 



a(t) 



4a (6a3ff 2 - 2 K p f /3 



2/3 



4(3H a (t - t ) + 4) 2 - 3k Pq (* - to) , (137) 



where we have used the initial conditions a (to) = ao, and 
a (to) = aoHo, where Hq = H (to). 

The general solution of the cosmological field equations 
(|128j) and (|129|) can be obtained as 



t 



1 



Aa^ 2 da 



B 



A- 10 da (138) 



^16nGp J V a o - B ' 

where ao is an integration constant and we have denoted 

A= (3ap +4a 3 ) 1/12 . (139) 

One can then obtain the Lagrange multiplier A in terms 
of a as 



A 



1 



2p a 3 



iapoH 2 + 6H 2 a 2 - 2np dt. (140) 



In classical mechanics, the Lagrange multiplier has the 
meaning of the force that keeps the constraint on the 
mechanical system. In the models with energy conserva- 
tion, we would like to conserve the energy of the ordi- 
nary matter, which amounts to provide some energy to 
the gravitational system. Eq. (|140[) can be written in a 
differential form as 



o^ Aa3 ) = 2^( 



P. 2 K P° 
6(1 - 2 lP 



(141) 



showing that the time variation of the Lagrange multi- 
plier density is proportional to the square of the Hubble 
parameter. 
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VIII. DISCUSSIONS AND FINAL REMARKS 

In this paper we have extended the work initiated in 
[301 ] and [3l| by considering a more general gravitational 
action in which the Lagrangian of the field explicitly de- 
pends not only on R and T, but also on the contraction of 
the matter energy-momentum tensor with the Ricci ten- 
sor. The gravitational field equations have been obtained 
in the metric formalism in two cases, corresponding to a 
non-conservative and conservative physical system, re- 
spectively. In order to impose the condition of the con- 
servation of the energy-momentum tensor we have used a 
Lagrange multiplier method, which implies the introduc- 
tion of a new vector field in the gravitational action. The 
equation of motion of massive test particles was derived 
in the non-conservative case and so was its Newtonian 
limit, corresponding to weak gravitational fields and low 
velocities. A density-dependent supplementary accelera- 
tion, acting on massive test particles, is induced in the 
presence of a non-minimal coupling between geometry 
and matter. The extra force on massive particles gener- 
ated by the geometry-matter coupling is always present, 
even in the case L m = p, and causes a deviation from 
geodesic paths. The presence of the extra force could ex- 
plain the properties of the galactic rotation curves with- 
out resorting to the dark matter hypothesis. It is inter- 
esting to note that this supplementary acceleration is also 
proportional to the matter density gradient, tending to 
zero for constant density self-gravitating systems. A sim- 
ilar dependence on the gradient of the Newtonian gravi- 
tational potential also appears in the generalized Poisson 
equation. 

The viability of the theory was studied by examining 
the stability of the theory with respect to local perturba- 
tions. In pure f(R) gravity, a fatal instability develops 
on time scales of the order of 10 26 s when the function 
f(R) satisfies the condition f"(R) < 0. This instabil- 
ity, called the Dolgov-Kawasaki instability, was discov- 
ered in the prototype model f(R) = R — [J. 4 /R, with 
(i ~ H ~ 10 ~ 33 eV [3g], which is therefore ruled out. 
In the present case, the condition of the stability with 
respect to the local perturbations can be formulated as 
!rr (Ro) - (po - To/2) f RT M. (R ) > 0, where i? is the 
background Ricci scalar. 

The cosmological implications of the theory were also 
investigated for both conservative and non-conservative 
theories. For this study we have adopted four functional 
forms for / (R, T, R ltv T liV ). In the non-conservative case 
we have shown that for two choices of the function /, 
the gravitational field equations admit an exponential, de 
Sitter type solution. Therefore matter-geometry coupling 
may be responsible for the late time acceleration of the 



Universe, as suggested by the observation of the high 
redshift supernovae [5|. An interesting solution of the 
field equations was obtained in the case of a conservative 
model with f (R,T, R^T^) = R + aR^T^. In this 
case if the coupling constant a > 0, the solution has 
an oscillatory behavior, with alternating expanding and 
collapsing phases. For a < 0, the scale factor of the 
Universe has a hyperbolic cosine type dependence. We 
have also investigated models containing the square root 
of the trace of the energy-momentum tensor. In this 
case in the high density limit the Universe has a super- 
accelerated expansion, but no de Sitter type phase can 
be obtained analytically. 

Work along similar lines has been done independently 
in [37[, although in a different setting, and with a dif- 
ferent focus, with mainly the cosmological aspects of 
the theory being investigated. Indeed, the authors of 
[371 ] mainly considered the accelerating solutions of the 
f(R, T, R^vT 11 ") theory, and attempted to find the func- 
tional form of / analytically. On the other hand we dealt 
with the other aspects of the theory, including the mo- 
tion of a test body in the gravitational field, as well as 
the Newtonian limit and the generalized Poisson equa- 
tion. We have also considered some cosmological solu- 
tions for the model. An important new result in our 
work is the use of the Lagrange multiplier method to 
implement en ergy - momentum conservation. In both our 
work and in [371 ]. the Dolgov-Kawasaki instability was 
explored, and the same results were obtained. 

The field equations of / (R, T, R liv T tiV ') gravity are ex- 
tremely complex. For different choices of the function 
/, cosmological solutions with many types of qualitative 
behaviors can be obtained. These models can be used 
to explain the late acceleration of the Universe, without 
resorting to the cosmological constant, or to the dark en- 
ergy On the other hand, this theory can open a new 
perspective on the very early stages of the evolution of 
the Universe, and may provide an alternative to the infla- 
tionary paradigm, which is facing very serious challenges 
due to the recently released Planck results. 
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